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The growth of structure in the universe begins at the time of radiation-matter equality, which 
corresponds to energy scales of ~ OAeV. All tracers of dark matter evolution are expected to be 
sensitive to neutrino masses on this and smaller scales. Here we explore the possibility of using 
cluster number counts and power spectrum obtained from ongoing SZ surveys to constrain neutrino 
masses. Specifically, we forecast the capability of ongoing measurements with the PLANCK satellite 
and the ground-based SPT experiment, as well as measurements with the proposed EPIC satellite, 
to set interesting bounds on neutrino masses from their respective SZ surveys. We also consider an 
ACT-like CMB experiment that covers only a few hundred deg 2 also to explore the tradeoff between 
the survey area and sensitivity and what effect this may have on inferred neutrino masses. We 
find that for such an experiment a shallow survey is preferable over a deep and low-noise scanning 
scheme. The precision with which the total neutrino mass can be determined from SZ number 
counts and power spectrum is limited mostly by uncertainties in the basic cosmological parameters, 
the mass function of clusters, and their mean gas mass fraction; all these are explicitly accounted 
for in our statistical Fisher-matrix treatment. We find that projected results from the PLANCK 
SZ survey can, in principle, be used to determine the total neutrino mass with a (la) uncertainty 
of 0.28eV , if the detection limit of a cluster is set at the 5a significance level. This is twice as large 
as the limits expected from PLANCK CMB lensing measurements. The corresponding limits from 
the SPT and EPIC surveys are ~ 0A4eV and ~ 0.12eV, respectively. Mapping an area of 200 deg 2 , 
ACT measurements are predicted to attain a la uncertainty of 0.61 eV; expanding the observed 
area to 4,000 deg 2 will decrease the uncertainty to 0.36 eV. 

PACS numbers: 



I. INTRODUCTION 

The temperature anisotropy and polarization state of the CMB constitute a unique window on the physical state 
of the universe in the era of recombination (z w 1100), ~ 400, 000 years after the big bang. Measurements with the 
COBE and WMAP satellites, along with a host ground-based and balloon-borne experiments, have led to impressively 
precise values of the basic cosmological parameters. Other cosmological probes, sensitive to certain combinations of 
these parameters, supplement CMB data, significantly reduce parameter degeneracies, and provide some control of 
systematic uncertainties. A major advantage of the CMB is its well understood underlying physics, a fact that 
justifies its distinction as the most precise cosmological probe. This said, it is also well recognized that the diagnostic 
power of the primary CMB anisotropy, whose main features reflect conditions at recombination, is rather limited 
in constraining certain cosmological parameters. Measurements with the recently launched PLANCK satellite arc 
cosmic-variance-limited up to i ~ 2500, and are therefore expected to reduce this limitation significantly by mapping 
the primary CMB anisotropy down to scales of a few arcminutes. 

Neutrino masses are a good example for the need to sensitively map the primary CMB anisotropy at high I. Below 
their Jeans scale, massive neutrinos free-stream and act as a hot dark matter (HDM) component, reducing the growth 
of the large scale structure (LSS) as compared with the evolution of structure in a pure cold dark matter (CDM) 
model. This characteristic neutrino signature provides a diagnostic measure that can be exploited to either determine 
or constrain their masses. Current cosmological constraints on neutrino masses come from the CMB temperature 
anisotropy e.g. [1,2], weak gravitational lensing [3-6], as well as galaxy [7-11] and Lya surveys [12-15]. WMAP7 data 
already set an upper limit (at the 95% confidence level) of 0.58 eV on total neutrino mass [2]. This limit, which 
comes mainly from CMB constraints on the early integrated Sachs Wolfe (ISW) effect (I < 200), is related to the 
fact that the neutrino temperature at radiation-matter equality is ~ 0.4 eV (in equivalent energy units), and that 
there are three active neutrino species. The ISW effect is associated with the kinematics of neutrinos at the time 
period between radiation-matter equality and recombination. This is an indirect probe of neutrino kinematics at, and 
prior to, last scattering. To improve on this bound, probes of neutrino kinematics at lower redshifts are required, 
where neutrino temperatures are lower (it can certainly be the case that a massive neutrino which is nonrelativistic 
today, was actually relativistic, i.e. freely streamed, at much higher redshifts). This may be achieved, for example, 
by employing lensing extraction of the CMB, e.g. [16-18]. 

The relatively strong bounds deduced from CMB data on the total neutrino mass should be contrasted with the 
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mass that would have provided closure density (of a hypothetical neutrino-dominated) universe, w 45 eV, and to the 
actually measured value of the matter density parameter, £l c d m ~ 0.3 (equivalently, ~ 5eV per neutrino species if the 
entire dark matter was made up of neutrinos) . From this perspective, the total neutrino mass is bounded to be lower 
than 4% (2a) of all dark matter (DM). As a tracer of LSS, gravitational lensing of the CMB is a sensitive probe of 
any cosmological parameters that affect the formation or growth-rate of the perturbed background matter density. 

Galaxy clusters are excellent probes of structure formation, and therefore can be used, in principle, for indirectly 
setting limits on neutrino masses through their dynamics and impact on structure formation. For this purpose, it 
has been suggested [19] that cluster number counts and cluster correlations can be used to set useful upper bounds 
on neutrino masses. These probes are sensitive to the evolution of structure down to the typical correlation scale of 
galaxy clusters, a few tens of Mpc, and are therefore sensitive to sub-eV neutrino masses. Lensing of the CMB with 
planned good sensitivity and high resolution experiments, e.g. [20,21] could do even better, perhaps even constrain 
the total neutrino mass, M v , at the 0.05 eV level. It is interesting to note that this level overlaps with the lower mass 
limit from neutrino oscillations, e.g. [22-24]. This implies that, to the level of confidence wc have in these forecasts, 
in the not too distant future, neutrino masses and mass hierarchy will be resolved, if next generation CMB satellites 
(e.g. EPIC [20]) will indeed be launched. However, other probes are certainly of interest even if their respective 
precision levels cannot compete with the CMB lensing probe, since they can provide much needed consistency checks 
for systematics and parameter degeneracy breaking. 

We have explored the possibility of determining neutrino masses from survey measurements of the SZ effect (Comp- 
tonization of the CMB by hot gas in clusters, e.g. [25-27]) by the currently operational ground-based SPT project, 
PLANCK satellite, and the proposed EPIC satellite mission. We also explored various scanning schemes for the 
ground-based ACT to gain some insight on how to optimize SZ surveys for neutrino mass inference. The basic moti- 
vation for doing so is provided by the fact that the SZ effect is a sensitive tracer of the evolution of the LSS. However, 
in contrast to the linear physics characterizing much of the evolution of the primary CMB anisotropy, the SZ effect is 
more model-dependent (e.g., mass-temperature scaling relation, gas mass fraction and spatial profile). A substantial 
source of uncertainty is the distribution of the entire cluster population in mass-rcdshift space, i.e. the mass function. 
This important function, whose specific shape and normalization reflect the details of the growth of density fluctua- 
tions, and the non- linear hierarchical collapse and mergers of sub-structures, can be best studied by state-of-the-art 
large-volume hydrodynamical cosmological simulations. Although very advanced, currently available computer codes 
predict a range of mass functions. At present, this mass function indeterminacy largely sets the precision limit of 
forecasting the total neutrino mass from cluster SZ number counts and power spectrum. 

In this paper we discuss constraints on neutrino masses that can be obtained from SZ surveys alone. In particular, 
we determine how likely it will be for PLANCK, SPT and EPIC to constrain neutrino masses by using the sensitivity 
of cluster number counts and SZ power spectrum to neutrino masses. For the first time in this context, attention is 
given to the robustness of the results to the assumed fiducial neutrino mass. One of our main objectives is to check 
to what extent the upper limits on neutrino masses derived from SZ surveys are robust to modeling uncertainties and 
other systematics. 

The paper is organized as follows. In section II we discuss the impact of massive neutrinos on the LSS with 
special emphasis on its manifestation in the transfer function and cluster mass function. Cluster number counts, 
correlation, and induced CMB anisotropy are described in section III. We lay out the methodology of our estimates 
for expected number counts and SZ power as a function of neutrino mass and address issues such as detection threshold, 
foregrounds, etc. along with a description of the Fisher matrix analysis used in this work in section IV. Our main 
results are presented in section V and further discussed in section VI. 

II. NEUTRINO IMPACT ON GROWTH OF THE LSS 

At a redshift 3200, the CMB and matter energy densities, together with the neutrino contribution, set the time 
of matter-radiation equality, t eq . This transition preceded recombination, but not by much; thus, the universe was 
not purely matter-dominated at recombination. Consequently, gravitational potential wells could decay and induce 
the early ISW effect, which boosts the angular power spectrum of the CMB temperature anisotropy on multipolcs 
that correspond to the horizon scale at recombination, i.e. I <~ 100 — 200. In other words, a rec /a eq , the ratio of the 
scale factors at recombination and matter-radiation equality is a measure of the linear growth of structure during the 
period between radiation-matter equality and recombination. This ratio clearly depends on neutrino masses. The 
fraction of total energy in the form of radiation (which would include neutrinos, if the total neutrino mass is below 
<~ leV) at recombination determines the amplitude of the early ISW effect. 

More stringent constraints on neutrino masses could, in principle, be obtained from the SZ effect. A massive 
neutrino that is nonrelativistic today [m„c 2 > kT v (0)], was relativistic at higher redshifts [m„c 2 < kT v (z)], and 
therefore could not easily cluster then, causing a relative suppression of structure formation compared to what would 
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be expected had this species been non-relativistic for the entire history of structure formation (t > t eq ). In principle, 
precise measurements of the LSS - for example, cluster number counts - can be used to place constraints on neutrino 
masses. 

The evolution of structure in the matter-dominated era is described in terms of the matter power spectrum, 

P m {k,z) = Ak n T 2 {k,z), (1) 

with the primordial density fluctuation spectrum, Ak n , where A is an overall normalization, n is the tilt of the power 
spectrum, and the transfer function, T(k,z) — T(M v ;k, z). An important quantity gauging the amplitude of the 
processed power spectrum observed today is the mass variance parameter on a scale of 8Mpc /i -1 , 

f°° Ah- 

a 8 = J P m (k,z)W 2 (kR)k 2 -^ 7 (2) 

where W(kR) is a window function, and R = 8Mpc h . The latter quantity incorporates the physics of neutrino 
damping; thus, a$ is a function of not only A, n, but also of neutrino masses (via the transfer function), and indeed 
any other cosmological parameter which affects structure evolution on scales of few tens Mpc and lower. Because these 
scales are comparable to typical scales of diffusion damping of density fluctuations by neutrinos with small masses, 
we expect that M v and u% will be anti-correlated. Since the SZ signature is a strong function of ag, it is expected to 
be sensitive also to M v . 



A. The Transfer Function 



Quantitative description of the effects of neutrinos on the evolution of the LSS begins with a modification of the 
density fluctuation spectrum as described by the transfer function. Whereas the standard CDM transfer function is 
separable in the perturbation wave-number k and redshift z 

T(k,z) = D{z)f{k) (3) 

where D(z) is the linear growth factor and, e.g. [28], 

f(k) = ln(1 2 + 3 4 g 34g) [1 + 3-89? + (16.1q) 2 + (5.46g) 3 + (6.7I<z)T 1/4 (4) 

and q = k/(fl m h 2 )Mpc^ 1 , a mixed DM model generally cannot be cast in such a simple form. This is due to the 
rcdshift-dcpcndence of the free-streaming scale, kFS, which is determined by both neutrino mass and the redshift- 
dependent neutrino temperature, T u (z). A reasonable first approximation to the transfer function in the presence of 
neutrinos, is e.g., [29] 

T 2 {k,z) -> D 2 {z)f 2 {k); k<k nr 

T 2 (k,z) D(*) 2 -*'"(l-8/„)T 2 (fc); k > k nr . (5) 

In the following this approximation is merely used in justifying order of magnitude sensitivity of the SZ power 
spectrum to changing the fiducial neutrino mass; in our numerical calculations we employ a publically available code 
[30] to precisely determine T(k,z). The effect of including nonlinear terms in f v on the matter power spectrum was 
investigated in several recent works, e.g. [31-33]. Indeed, this higher order effect modifies the matter power spectrum 
on scales of few Mpc and smaller, and in principle should be accounted for when data from future SZ surveys is 
available. In the current work this effect is ignored for simplicity. We can justify doing so by noting that the extra 
0(f 2 ) term will only increase the sensitivity of a 8 , and consequently also of the SZ effect, to the neutrino mass and 
thereby will increase the sensitivity of SZ-based cluster number counts and power spectrum to neutrino mass. In 
this regard the present analysis represents a conservative assessment of the capacity of cluster number counts and 
cluster correlations to constrain neutrino mass. Indeed, the latter studies [32,33] show that the factor 1 — 8/„ should 
be replaced by « (I — 9.8/„) for k 3> k nr . It is also clear that the f 2 corrections to the matter power spectra are 
important when /„ is sufficiently large. Indeed, neutrino masses twice or three times as large as the current 95% 
confidence upper limit on neutrino mass from WMAP (0.58 eV) were considered [31-33]. For the range of neutrino 
masses considered in this work, M v < 0.6 eV, the nonlinear neutrino effects are smaller. 
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B. Cluster Mass Function 



The basic quantity which describes the number density of clusters as a function of their mass and redshift - basic 
properties by which clusters are identified - is the mass function. As will be shown later, the total neutrino mass can 
be derived from comparison of the observed number of clusters at a given redshift-bin to the number predicted from 
the mass function, d "^ ;z - ) . The latter function is defined in terms of the differential cluster number in a redshift 
interval, 

dN{M, z) = f sky M ^ z) dVdM, (6) 

where dV is a volume element and f s k v is the observed sky fraction. The total number in a given interval Az around 

Zi is 

A* (*) = U y Az^ J **^dM. (7) 

As noted earlier, the currently most optimal determination of the mass function is from extensive cosmological 
simulations. Here we adopt the mass function derived recently from a large set of collisionless cosmological simulations 
in ACDM cosmological model [34] . Expressing the mass function in the familiar form, 

dn p m d\n.{a~ 1 ) 



these authors derived a fit of the form 



(9) 



where the parameters A, a, b and c were best-fit to the results of their simulations. The mass function does not have 
a universal form, a fact that becomes apparent by the deduced dependence of these fit parameters on both redshift 
and the overdensity at virialization, A v . 



III. SZ POWER SPECTRUM 



The SZ effect is a unique probe of cosmological parameters and cluster properties; its statistical diagnostic value 
is reflected through number counts and the power spectrum of the CMB anisotropy it induces. The dependence of 
cluster number counts and cluster correlations on the total neutrino mass, and an assessment of the feasibility of 
actually determining it from their measurements, are our main objectives in this work. To this end, we carry out a 
precise calculation of the SZ signature imprinted by a cluster of mass M at redshift z. We then calculate the integrated 
statistical pattern of the SZ-induced CMB anisotropy projected on the sky, and the corresponding cluster number 
counts. Obviously, the ultimate practical issue is how well can clusters be detected by a specific CMB experiment 
(characterized by both detector noise and angular resolution) in the presence of the primary CMB anisotropy and the 
emission of foreground (point) sources. 

The thermal component of the SZ effect induced by a cluster with optical depth r = / n e OTdl and gas temperature 
T e induces a small level of fractional change in the CMB temperature [35,36] 

AT 

-Y = vH*) (io) 

where y = rO e is the Comptonization parameter, e = kT e /(m e c 2 ), and k, m e , and c are the Boltzmann constant, 
electron mass and speed of light, respectively. n e and <jt are the gas number density and Thomson cross section, 
respectively, and the integration over the photon pathlcngth is along the line of sight. In the classical limit the spectral 
function reads 



T{x) =xcoth(x/2)-4 (II) 

where x = hu/(kT) is the dimensionlcss CMB frequency, and h and T denote the Planck constant and CMB temper- 
ature, respectively. Note that in this work we ignore relativistic corrections [37, 38], whose weighted impact over the 
full temperature range of clusters is small. Since n e is non-uniform and has a typical scale, the cluster core radius r c 
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(as well as its temperature T e which we here assume uniform T e (r) = Tq for simplicity), the projected 2D profile of 
the y-parameter on the sky, inherits this typical angular scale, that determines together with other parameters (as we 
discuss below) the SZ flux associated with the cluster and thereby the likelihood of its detection with a given CMB 
survey. The 2D Fourier transform of the projected 3D Comptonization parameter generated by an individual cluster 
is, e.g. [39] 
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where 



(1a/ V c and the angular diameter distance in a flat ACDM model is 

dz' 

d A {z) = 



{l + z)H J y/n m (l + z'^ + SIa 



(12) 



(13) 



with H , fl m and denoting the Hubble constant today and matter and vacuum energies in critical density units, 
respectively. Note that the neutrino energy density is contained in £l m since all neutrino masses considered in this 
work are certainly non-relativistic out to the highest- redshift clusters considered in this work, z ~ 1. The virial 
radius of the cluster is related to its virial mass via the relation M v = ^-R^A v {z), where A v (z) is the overdensity at 
virialization with respect to the background, calculated from the spherical collapse model in ACDM universe. 

An expression for the virial radius, R v , as a function of the virial mass, M v , and the underlying cosmology can be 
derived noting that 



M v = ^Rl Pv {z) = ^RlA v (z)^ )Pc (0), 



(14) 



where p c (z) and A c (z) = p v {z) j p c (z) denote the critical density and the virial overdensity at redshift z, respectively. 
This relation leads to 
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where Mi 5 is the cluster mass in units of 10 M©. Writing 

Pc(0) 



Hi 



and the explicit redshift dependence of the Hubble parameter 

h 2 (z) = h% [n A + n m (i 

we obtain 

gc(0) _ 1 n m (i+zf 

pi z \ n A +n m (i+z) a fi A +f2 m (i+z) 3 f2 m (i+z) a — 

Substituting this relation in the expression for the virial radius leads to 
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Hot intracluster (IC) gas is assumed to follow a (3 density profile, 

n e (r)=n e , [l + (r/r c ) 2 ]- 3 ^ 2 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



with the observationally-deduced j3 — 2/3, and a core radius defined as 10r c = R v . The central electron gas density 
is n e fl = Pgfl{X + l)/(2m p ), where X and m p denote the Hydrogen mass fraction (taken to be 0.69) and the proton 
mass, whereas 



Pg,o 



fg-M v 



47rr3/ 10 ' x 2 {l + x 2 )-W 2 dx 



(21) 
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is the central gas density, assuming that IC gas constitutes a fraction f g = 0.12 of the total virial mass. The virial 
expression is adopted for the gas temperature 

k B T e , -— — , (22) 

valid for R v r c , where /i = 0.59 is the mean molecular weight. The resulting dependences of the core radius, gas 
density and temperature on the cluster mass and redshift obviously stem from these definitions. 

A key quantity employed in the next section in the calculation of the signal-to-noise detection of an individual 
cluster of mass M at redshift z is the power spectrum generated by the cluster 

Cf z {M,z-M v ) = ±-\£(l,M„;M,z)\ 2 T 2 (x). (23) 

Clearly, the total contribution from the entire cluster population to the SZ signal (which is irrelevant to this work) is 
obtained from Eq.(23) by a weighted integration over M and z with ^j^z"' 1 - e - integration over differential cluster 
number dN (Eq.6). The SZ power spectrum is shown in Fig.l for the three fiducial neutrino masses considered in this 
work: 0.1 (black), 0.3 (blue), and 0.6 (red) eV. All power spectra were obtained by setting T = — 2 which corresponds 
to the Ray leigh- Jeans limit. The appreciable sensitivity to the assumed neutrino mass comes from both the strong 
dependence of the transfer function on /„ (Eqs.2, 1 & 5), i.e. Aas/as ~ —8Af v and the scaling of the SZ power 
spectrum with high power of as, typically a\. 

The contribution to the SZ power spectrum from cluster correlations is typically an order of magnitude lower 
than the Poissonian component, and peaks on angular degree scales where it is overwhelmed by the primary CMB 
anisotropy. Consequently, its expected relative contribution to the angular power spectrum is small (e.g., [40, 41]) 
and we therefore do not consider it in our analysis. Also ignored is the small kinematic component of the SZ effect 
[36]. Had we included these two components in our analysis, cluster detection rate would have increased somewhat, 
improving our statistics and therefore lowering the uncertainty on neutrino masses, but these effects are very small. 

IV. FISHER MATRIX ANALYSIS 

Our cosmological model includes the normalization A and tilt n of the primordial scalar perturbations, neutrino, 
fl^/i 2 , n m h 2 , and baryon, flbh 2 , energy densities, the Hubble parameter h in km/sec/Mpc units, dark energy equation 
of state w, optical depth for reionization r, and the primordial helium abundance Y p . The last two parameters do 
not directly affect the SZ quantities but are correlated with other cosmological parameters which might affect the 
uncertainty of these and indirectly impact neutrino mass when probed with SZ number counts (and power spectrum). 
We first carry out the forecasts for the primary CMB. This is done both with and without lensing extraction (LE). 
The corresponding Fisher matrices are denoted F pr and F LE , respectively. 

The extremely steep dependence of the SZ effect on as is clearly expected to result in strong dependence of neutrino 
mass constraints on the fiducial neutrino mass. This follows from the fact that for lower neutrino masses (when A is 
fixed) as increases and as a result the SZ amplitude increases, boosting the signal-to- noise (S/J\f) of cluster detection 
and tightening the constraints on inferred neutrino mass. Clearly, the main cluster parameter which determines the 
magnitude of the SZ signal, and thereby whether or not it will be detected in a flux-limited survey, is the total mass. 
However, the cluster redshift is also of relevance, as explained above and will be further discussed below. Since the 
mean properties of the cluster population provide a good basis for estimating global parameters, their gaseous mass 
fraction could, in principle, be related directly to the ratio Q;,/Q TO (after accounting for the small contribution by the 
cluster member galaxies). However, in this work we adopt the fiducial value of the gas mass fraction deduced from the 
sample of SZ measurements with the BIMA & OVRO telescopes, f g = 0.12 [42]. Due to the appreciable uncertainty 
in the measured value of f g this cluster parameter is also treated as a model parameter in our Fisher matrix analysis. 

In calculating the Fisher matrix for the primary CMB (with and without LE) we follow the standard approach 
which we do not reproduce here; details of the calculation can be found in (e.g.) [17]. 

A. M v from Cluster Number Counts and Correlations 

As in, e.g., [43], we write for the Poissonian likelihood function for cluster number counts in the i'th redshift bin 



(24) 
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where Nf is the predicted (for a given cosmological model and specified neutrino mass) and N° the observed cluster 
number in the i'th redshift bin. Note that here and in the following we use abbreviated notation and denote AiV, 
of Eq.(7) as Ni. If Nf is a function of several cosmological parameters Afc, we may consider small deviations with 

respect to the expected fiducial value N° w Nf + ^2 k ^or~AAfc. The Fisher matrix for cluster number counts can then 
be determined from the likelihood function 

N = a 2 \nc = v j_mm . , 

Jm d\jd\ m Y N < dX i dX ™ 
The estimated uncertainty in the parameter Aj is then 

M i = (F? j )- 1/ *, (26) 

where we take the square root of the j'th Fisher matrix element. 

The second neutrino probe we employ in this work are measurements of the SZ power spectrum. Again, only those 
clusters which are detected at 5a or better qualify for our 'simulated' sample. For the corresponding Fisher matrix 
we use the Fourier transform P c (k; z) of the cluster correlation function C(r; z), as in [44] 

Pc = v dHklk\\Pc)ai dHkiWai (VkVeff(k))ai ( 

jm d\j d\ m ■> 1 ' 



a.i 



where the effective survey volume is 



/dr 
r2<yZ ^Tz dz 



n( Zi )P c (k) 1 



1 + n{ Zi )P c {k) 

and n{zi) is the number density of detected clusters at the i'th redshift bin. V4 is the k-space volume element 



(28) 
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P c {k^, fc||; z), the cluster power spectrum, is 

P c (k ±7 k ll ;z) = [l + f3(z)^] 2 b 2 (z)P(k;z) (30) 

where k 2 = k\ + k 2 , \x = k\\/k, fi{z) = d ^^a)^ an< ^ ®9 r ^ s ^ e nnear growth factor obtained from solving the 
equation governing the linear evolution of density perturbations 8 

d 2 8 \2 1 /ft m a- 3 -2Q A \l 8 3 Q m (l - f v ) 

-o = (J (31) 



da 2 



a 2a \ ft m a 3 + Q\ 



da 2 fl m a 2 + il\a 5 



where a is the scale factor, and /„ = fl v /Q m is the fraction of DM in the form of neutrinos. Here, the mass-averaged 
linear bias is 



b(z) 



where b(M; z) is obtained from hydrodynamical simulations [45]. 

The total Fisher matrix that includes cluster number counts combined with either the primordial or lensed CMB is 
Fjm + Fjln or Fjm + ^jli respectively, and the uncertainty in the parameter Xj is similar to that in Eq.(26), except 
for the fact that now is replaced with either Fj^ + Fj^ or Fj^ + Fj^ . Similarly, from the power spectrum we 
obtain Fj^ + F? m or F^ + F? m , respectively. 

Integrated cluster quantities, such as number counts and the power spectrum, obviously depend on the assumed 
mass range. Our nominal mass range is taken to be 10 13 — 3 x 10 15 M . While it is doubtful whether our description 
of IC gas properties is valid at the low mass end of this range, we have verified that our results are insensitive to the 
exact value of the low mass end by repeating the calculations with three different values for the low mass cutoff, with 
only a negligible fraction ( 2% and 0.2% for EPIC and SPT, respectively) of the clusters in our mass range have 
masses smaller than 5 x 10 13 ft. _1 M Q . This is an important consistency check since the assumption of virialization 
(Eqs. 14-22) is expected to break down at around this value, and since the cluster mass cannot be deduced in a 
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model-independent way, it is reassuring that the detection threshold is set such that clusters for which our model 
breaks down are not included in the sample. This essentially is the case as long as the detection threshold is taken 
to be at the 5a significance level. The upper mass limit of this range is in accord with recent detailed joint data 
analyses (e.g. [46]) that have resulted in reliable mass estimates that are already very close to the above value. Here 
too, raising the upper mass to 10 16 M Q has little impact on the uncertainty in inferred M v . 

To estimate the S/Af with which a cluster can be detected in a survey, we assume that the main noise sources 
are the instrument, primary CMB anisotropy, and point source emission. We ignore other astrophysical foregrounds, 
such as free-free emission, synchrotron radiation and dust, due to their Galactic origin and the planned removal of 
this sky region in the PLANCK and EPIC data analysis. Galactic contamination in the SPT and ACT survey is 
minimized by selecting a radio-quiet sky region. Optimal (minimum variance) unbiased filters can be constructed for 
cluster detection with the resulting S/Af variance [47, 48] characterized by minimum confusion noise 
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which is merely proportional to the ratio of the power spectra of signal and noise weighted by Fourier space volume 
(oc I) and summed over all accessible multipolcs. Contribution to the S/Af from scales below the beamsize are 
exponentially damped 

Cf et = (A T 6 )y 2<7 > (34) 



where 6 b = ^/8rn(2)er&. At is a noise-equivalent-temperature (NET) of the experiment, and the sum over i runs over 
all the detector frequency channels. The relevant parameters for the PLANCK instrument, SPT, and EPIC are listed 
in Table I. The adopted power spectrum of point sources is [49] 

C[ s = A ps (x Q /x)- 4 (e x - l) 4 x- 4 e- 2 *, (35) 

where A ps = 0.015^ 2 -str, x Q = (hv Q )/{kT), and v Q = AQ.lGHz is WMAP's Q-band. Cf s , A T , and b are all 
frequency-dependent. 

As shown in Eqs.(12) & (23) the SZ power spectrum is a function of the cluster core radius, 8 C , which in itself is a 
function of M and z (Eqs. 19 & 20). Therefore, one can define the S/Af per a single cluster in the M-z cell 

S\ 2 _^ fldl \it{M,z;vJ\ 2 

(36) 



m, z " Y J 2lT W et ^ + c f MB + °f 3 ^ ' 

where £;(M, z; Vj) is given by Eq.(12); recall though that the cluster core radius is a function of both mass and redshift. 
This measure, together with the mass function, which determines the number density in each of the M-z cells, can 
be used to estimate how many clusters will be detected at a given cell by a specific experiment. This criterion is 
employed in the next section. 



V. RESULTS 



In this work we have adopted the ACDM cosmological model with WMAP best-fit parameters. The cluster popula- 
tion is described in terms of the mass function of Tinker et al. [34, 45], and the (hot) IC gas mass fraction was taken 
at the level deduced from SZ observations, f g ~ 0.12. IC gas has been assumed to be isothermal with a beta density 
profile. It is important to note that in this analysis no priors were set neither on the cosmological parameters nor on 
f g . We reiterate that changing neutrino masses over the range [0.1, 0.6] eV changes as by ~ 25%. Since the SZ effect 
strongly depends on as, we expect observational measures of the effect to show a corresponding level of sensitivity to 
varying M v , as is indeed evident from Figure 1. 

Our analysis includes several inherent noise sources - detector noise, and the noise due to the primary CMB 
anisotropy and extragalactic point sources. We first calculate the distribution of (S/N)m,z with which clusters of 
mass M at redshift z will be detected (Eq. 36) in the PLANCK, SPT and proposed EPIC surveys. As is clear 
from Eq.(36) the clusters detectability likelihood depends on both their masses and redshifts; the mass and redshift 
together determine the S/Af, as is clear from the fact that the gas temperature and core radius (both physical and 
certainly its angular projection on the sky) depend on both M and z (Eqs. 19 & 20). We set the level of statistical 
detectability to 5a. 

Our basic results for neutrino mass uncertainty based on cluster number counts and power spectrum from SZ 
surveys with the ongoing PLANCK, SPT, and proposed EPIC projects are presented in Tables II and III for the 
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fiducial neutrino masses, 0.1, 0.3, and 0.6 eV. In Table II we show (from left to right) the expected uncertainty on the 
inferred M„ from the primordial CMB (both temperature anisotropy and polarization), primary CMB and cluster SZ 
number counts, primary CMB and cluster correlation, and finally the constraints from joint primary CMB, cluster 
number counts and power spectrum. Table III is the same as Table II, but the priors are from LE. It is particularly 
interesting to see that if M v is ~ O.leF, then the joint primordial CMB, cluster abundance and correlation analysis 
for PLANCK yields a Mv = 0.28eV, whereas the constraint from CMB LE is O.UeV. 

Similar results for ACT are shown in Table V. This project is viewed here as a generic ground-based SZ survey 
experiment with deep exposure on a very small patch of the sky. Our goal is to optimize the scanning strategy for 
neutrino mass inference. The analysis presented here clearly illustrates that the uncertainty on neutrino masses can 
be reduced by ~ 1/2 when the survey area increases from 200deg 2 to 4,000deg 2 . Indeed, increasing f sky results in a 
higher noise per sky-pixel, which lowers the detection significance of many clusters, but with the 5a detection-threshold 
we adopt here it turns out that only the smallest clusters are removed from our sample and this is outweighted by 
improving the statistics of massive/nearby clusters by virtue of increasing f s ky Specifically, we considered four fiducial 
ACT modes; ACT-A, B, C and D (Table IV). ACT-A refers to the nominal ACT experiment that covers 200deg 2 of 
the sky. ACT-B, C and D scan 400, 1000 and 4,000 deg 2 , respectively. The larger f s k y , the shorter is the integration 
time available per sky-pixel, and assuming the integration is uniform over the sky and that the instrumental noise is 

— 1/2 

inversely proportional to the integration time, At oc t int , we can rescale ACT-A to the other three modes. While 
ACT-A has very low noise it is limited by the lower statistics because of the survey area, 0.5% of the sky. Results 
listed in Tables V and VI clearly illustrate that eventually the gain in sky area (in the cases of ACT-C and D) result 
in a notable improvement in neutrino mass constraints. The lesson from this is that observing a larger sky-patch is 
preferable, at least for neutrino mass inference and for a CMB telescope with ACT-like specifications. It should be 
noted here that the priors, i.e. constraints that come from either primordial CMB or LE that we adopted in these 
tables come from PLANCK, not ACT itself (which should also be affected by changing f s k y ), and are for a fixed 
fsky =0.65. 

An important aspect of the theoretical treatment is careful selection of the mass range. In principle, we could 
compare expected and observed number counts on the entire 2D M-z space, if it was realistic to reliably deduce 
cluster masses (from the 'blind' surveys considered here). Since cluster mass inference is highly model-dependent we 
always integrate the cluster distribution over mass and consider number counts in redshift-space only, as described 
earlier (Eqs. 25 & 27). Integrating the expected number counts and cluster correlation over mass we have to impose a 
low mass cutoff, thereby removing all low-mass clusters from our sample. As noted already, this is motivated mainly 
by the fact that objects that have detectable mass of sufficiently hot (> 1 keV) gas, and for which our theoretical 
description is still valid - including aspects such as the gas density profile, mass-temperature scaling, and indeed the 
attainment of hydrostatic equilibrium in the first place - do not include very poor clusters and groups of galaxies. 
We have checked and verified that our results are robust with respect to changing this low mass cutoff to the lowest 
cluster mass of 10 13 M© we consider in this work. This obviously reflects the consistency between our theoretically 
selected value of the low mass cutoff and the high significance that we require for the threshold for SZ detection. 

Formation of clusters is relatively late in the ACDM model, with the vast majority of medium to high mass clusters 
forming at z < 1 (as is clear from plots of the probability distributions of formation times in [46]). We checked that 
essentially irrespective of the fiducial neutrino mass considered, the contribution to the Fisher matrix saturates at 
Zmax ^ 1- Higher-redshift clusters are simply not sufficiently massive and in addition project on too small angular 
scales to be detected. 

Redshift binning may also introduce some bias; by increasing bin width the diagnostic capability is weakened, but 
the impact of random fluctuations e.g. modelling errors that bias the redshift of cluster formation is reduced. This 
motivated us to determine an optimal redshift-binning in our number counts analysis by binning the redshift interval 
[0 — 1] into 10, 20 and 40 linearly spaced redshift shells. We found that increasing this number from 10 to 20 improves 
the statistical error on the deduced neutrino mass by not more than ~ 30%. Further binning refinement to 40 shells 
reduces the error by at most 15%, implying that further increase in the number of shells does not reduce the statistical 
uncertainty but may come at a cost of a higher level of systematic error. All results reported in this work using cluster 
number counts were obtained by binning the simulated data into 40 redshift bins. For the cluster correlation probe 
we used only 10 redshift shells; as in [51] we choose k- binning in Eq.(27) to be wide enough for different bins to be 
uncorrelated. We assume Az = 0.1 for this analysis and Ak — O.OObhMpc^ 1 in the k-range [0. 005-0. V^hMpc^ 1 . 

The most important source of uncertainty in modeling cluster abundance and internal properties is the mass 
function. Current uncertainties in this basic function were explicitly included in our analysis. Continued extensive 
cosmological hydrodynamical simulations (e.g. [52]) are likely to result in a significantly more precise description of 
this function across the cluster mass range. In contrast, uncertainties stemming from using simple models for the 
spatial profiles of the gas density and temperature are of secondary importance, simply because these are much reduced 
when calculating integrated SZ measures and, more generally, will have little effect on cluster detection. After all, the 
magnitude of the effect in a cluster is not determined by each of these quantities separately (in the non-relativistic 
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limit that is valid for our purpose here), but rather by the integrated gas pressure along the line of sight. Thus, even 
though central regions of clusters are roughly isothermal, and the gas density profile is well approximated by the beta 
profile, realistic deviations from isothermal beta profile for the pressure, when integrated across a cluster and over the 
mass function, introduce a level of uncertainty that is ignorable for our purposes here. Mass function normalization 
and three additional shape parameters (Eq. 9) are included in our analysis. 

VI. DISCUSSION 

The CMB is perhaps one of best understood and model-independent underpinnings of modern cosmology. Temper- 
ature anisotropy and polarization measurements of the primary CMB signal already taught us a great deal about the 
geometry of the universe, its energy content, reionization history, etc. With the imminent detection of the lensing- 
induced signal of the CMB by PLANCK, and possibly even the inflationary- induced B-mode polarization (if inflation 
is indeed related to GUT scale physics), the CMB may open new windows into the physics of neutrinos as well as 
processes that occurred in the GUT era; the latter are well beyond the reach of any present or future terrestrial 
experiment. 

However, the capacity of primary CMB alone to constrain ~ Mpc cosmology is rather limited; even CMB lensing by 
the LSS, a sensitive probe of neutrino masses, takes place on considerably larger physical scales. In addition, standard 
forecasts of LE performance widely build on the assumption that the CMB signal is gaussian in the absence of lensing, 
and any nongaussianity present in the data can be attributed to lensing of the CMB. It should be emphasized that 
the results for neutrino masses derived from CMB and LE reported here, and elsewhere in the literature, always 
make this simplifying assumption. However, non-gaussianity can also be induced by astrophysical sources. The issue 
of primordial nongaussianity is still open and should be considered as another source of confusion. On the relevant 
angular scales secondary CMB signals, such as the SZ effect, are known to be present. These are typically nongaussian 
and can interfere with LE. 

Structure on Mpc scales probes the entire evolution history of matter perturbations down to these scales. This is 
especially relevant to neutrino physics via the effect of neutrino free streaming on these and larger scales. Indeed, 
Lya observations can supplement the CMB on these scales but since the astrophysics of the former is more involved 
it is desirable to consider other probes of clustering on Mpc scales to complement Lya, as well as other traditional 
probes such as galaxy clustering, shear measurements, etc., or at the very least to serve as a consistency check. 
The possibility of using cluster catalogues in constraining neutrino masses has already been discussed in [19]. More 
generally, extraction of cosmological parameters, such as fif,, ^a, from cluster number counts - in conjugation with 
other cosmological probes - is widely discussed in the literature, e.g. [43,53]. 

In this work we considered the prospect of SZ surveys to constrain neutrino masses. We have shown that using 
cluster number abundance alone neutrino mass uncertainties may be constrained to the ~ 0.36 — 0.55 eV levels by 
PLANCK and ^2 — 3 times weaker by SPT SZ surveys alone, on a wide range of allowed neutrino masses {M u < 0.6 
cV). The proposed EPIC satellite with its SZ survey of thousands of clusters will be able to set impressive constraints 
on neutrino masses at the level that will enable us to constrain models of inverted mass hierarchy (where M u > 0.1 eV). 
For this seemingly promising potential to be realized a full account of systematics has to be achieved; we elaborate on 
mass-function uncertainties below. Systematics, such as unknown gas fraction f g , are easily accounted for by simply 
including them as free parameters in our analysis; with predicted hundreds (SPT) or thousands (PLANCK and EPIC) 
clusters we have sufficient freedom to include other parameters, such as f g , in addition to M v and A. It is quite likely 
that follow-up X-ray observations will yield a more precisely determined value of f g , a possibility we clearly do not 
consider in this work as we assumed no prior on f g . 

Another subtlety comes from possible degeneracy with fide- Including degeneracy with dark energy in our framework 
is indeed possible with the code employed here for the transfer function [30] as it was developed especially for the 
purpose of calculating the impact of large scale structure suppression by both neutrino free streaming and accelerated 
expansion of the universe induced by dark energy. A very crude estimate of this effect can be obtained by fixing the 
CDM component, Sldm, and changing fl v on the expense of fide, but since the total energy density is assumed to be 
fixed to its critical value, and since both neutrino (with a total mass in the range considered here) and dark energy 
do not clump on galaxy cluster scales, we can conclude that - at least to first approximation - allowing all three 
fide and fldm to vary (subject to the constraint that the total energy density is fixed) does not substantially affect 
our results. 

The impact of sub-eV neutrino masses will obviously be reflected also in the superposed SZ power from the popula- 
tion of clusters. We have calculated the full power spectrum by integrating the expression for a single cluster (specified 
in Section III) over the mass function for three neutrino masses, 0.1, 0.3, and 0.6 eV. Results of this calculation are 
shown in Fig. 1 by the black, blue, and red curves, respectively. Since as is approximately proportional to 1 — 8/^, 
changing neutrino mass in the range [0.1 — 0.6] eV corresponds to varying /„ in the range [0.7 — 4.3] x 10~ 2 ; this 
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introduces a fractional change Aos/os ~ — 8A/„. The steep dependence of SZ power on erg translates to a relative 
change of « 56/„. From Fig. 1 we see that this factor amounts to an overall variation of the peak power by a factor of 
~ 3 for the above range of neutrino masses. Note though that in our calculations we used the precise transfer function 
and did not employ this crude approximation which we mention here only for the sake of obtaining a rough estimate 
of the SZ sensitivity to changing neutrino masses. Confusion with the CMB primary and other confusing signals 
renders this neutrino signature less diagnostically useful as compared with the more optimized cluster SZ number 
counts or correlations. It is important to note that the steep SZ power spectrum dependence on neutrino masses does 
not directly imply similar sensitivity in number counts. The only direct effect of the power spectrum dependence on 
M v is that large neutrino mass decreases the SZ power spectrum, and therefore the S/Af with which galaxy clusters 
are detected, thereby boosting the uncertainty on neutrino masses as is illustrated in Table II. 

Another tantalizing possibility is using these large cluster surveys to simultaneously calibrate the mass function and 
determine neutrino mass. This will help reducing possible bias that could be caused by miss-calibration of the mass 
function. In our analysis we adopted the functional form of the Tinker mass function, Eqs.(8) & (9), and considered 
the shape parameters as free parameters in the statistical analysis, adopting the central values reported in [34]. To 
test the robustness of our results to changing the mass function used in the parameter estimation we repeated the 
analysis with a different mass function which was similarly derived from hydrodynamical simulation - the Jenkins 
mass function [54]. We found that the constraints on neutrino mass do not change by more than 4%, and in most 
cases the differences do not exceed the 1-2% level. 

It is clear that more can be done to refine the analysis presented here since, admittedly, cluster physics and statistical 
properties of the cluster population are model-dependent. However, we believe that the adoption of a rigorous criterion 
for cluster detection in M-z space, assuming implementation of well-known unbiased optimal filters (in conjugation 
with the need to reject low mass clusters from our sample), the realization that the statistics improve when neutrino 
masses are lower as as increases, the use of a high precision transfer function for neutrinos in our theoretical SZ model, 
and optimal binning in redshift space, all these factors together constitute important steps towards achieving a high 
precision probe of neutrino masses. These features lend credibility to the results presented in this work. The results 
are encouraging: Cluster number counts, and to a lesser extent cluster correlations, may indeed be reliable neutrino 
probes if carefully analyzed. The utility of this probe largely depends on the neutrino mass; the lower the mass, the 
tighter is the constraint (as is shown in Table II and explained above). The prospects will only be further enhanced 
when SZ surveys will be jointly analyzed with results from extensive X-ray surveys. 
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FIG. 1: SZ power spectrum for three neutrino masses: 0.1 (black), 0.3 (blue) and 0.6 eV (red). As explained in the text, the 
dependence erg tx 1 — 8/„ translates to a relative change of ~ 56/„ in the value of as when the total neutrino mass assumes 
values in the range 0.1 — 0.6 eV. 
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Experiment 


/sky 


v\GHz\ 


HA 


A T [fiK] 






30 


33 


4.4 






44 


23 


6.5 






70 


14 


9.8 






100 


9.5 


6.8 


PLANCK 


0.65 


143 


7.1 


6.0 
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5.0 


13.1 






353 


5.0 


40.1 






545 


5.0 


401 






856 


5.0 


18300 






150 


1.0 


13.4 


SPT 


0.10 


274 


0.52 


71.4 






345 


0.29 


583.9 






30 


28 


0.50 






45 


19 


0.30 






70 


12 


0.21 






100 


8.4 


0.22 


EPIC (4K) 


0.65 


150 


5.6 


0.25 






220 


3.8 


0.66 






340 


2.5 


2.24 






500 


1.7 


9.41 






850 


1.0 


740.0 



TABLE I: Sensitivity parameters for Planck, SPT, and EPIC: / s k y is the observed fraction of the sky, vo is the central channel 
frequency (in GHz), Ob is the FWHM (Full- Width at Half-Maximum) in arcminutes, and At is the temperature sensitivity 
per pixel in /iK. Polarization sensitivity is Ae = As = \/2At, except for the 545 and 857 GHz bands of PLANCK which 
are unpolarized. We assume all the frequency bands can be used for primordial CMB and LE, but for the purpose of cluster 
detection we assumed that only the 100, 143 & 353 GHz (PLANCK) and 100, 150 & 340 GHz channels (EPIC) can be reliably 
used, since these are the least foreground-contaminated bands. 



Experiment 


M v [eV] 


<tm„ [eV](prim.) 




eV r ][prim.+N(z)] 


om v [eV] [prim.+Pc(z)] 


om v [eV][prim.+N(z)+Pc(z)] 




0.1 






0.46 


1.21 


0.44 


SPT 


0.3 


2.02 




1.25 


1.37 


1.01 




0.6 






1.52 


1.55 


1.27 




0.1 






0.36 


0.55 


0.28 


PLANCK 


0.3 


0.67 




0.55 


0.60 


0.42 




0.6 






0.55 


0.63 


0.50 




0.1 






0.13 


0.23 


0.12 


EPIC 


0.3 


0.34 




0.25 


0.26 


0.20 




0.6 






0.30 


0.30 


0.25 



TABLE II: Statistical uncertainty on total neutrino mass from cluster number counts obtained from PLANCK, SPT, and EPIC 
SZ surveys (no LE). Shown are neutrino mass uncertainties (lc) expected to be obtained from the respective planned surveys. 
We consider om v assuming fiducial masses 0.1, 0.3 and 0.6 eV. We show the results from CMB alone, CMB+cluster number 
counts N(z), CMB+cluster correlation Pc, and CMB+N(z)+Pc. 
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Experiment 


M v [eV] 


a Mv [eV]{LE) 




[eV][LE+N(z)] 


<tm„ [eV][LE+Pc(z)] 


om v [e^][LE+N(z)+Pc(z)] 




0.1 






0.35 


0.41 


0.25 


SPT 


0.3 


0.52 




0.50 


0.44 


0.34 




0.6 






0.51 


0.49 


0.43 




0.1 






0.12 


0.12 


0.11 


PLANCK 


0.3 


0.14 




0.12 


0.13 


0.12 




0.6 






0.13 


0.14 


0.13 




0.1 






0.04 


0.04 


0.04 


EPIC 


0.3 


0.05 




0.04 


0.04 


0.04 




0.6 






0.04 


0.04 


0.04 



TABLE III: Same as in Table II but with LE priors. 



Experiment 


deg 2 


v[GHz] 


[1'] 




ACT-A 


200 


145 


1.7 


2.0 






280 


0.93 


8.8 


ACT-B 


400 


145 


1.7 


2.83 






280 


0.93 


12.45 


ACT-C 


1,000 


145 


1.7 


4.47 






280 


0.93 


19.68 


ACT-D 


4,000 


145 


1.7 


8.94 






280 


0.93 


39.35 



TABLE IV: ACT-A, B, C and D: ACT-A is rescaled to B, C and D by assuming f sky oc tj„ t and A T oc t. 



Experiment 


o~m„ [eV](prim.) 


<?m„ [eV][prim.+N(z)] 


fM„ [eV] [prim.+Pc(z)] 


(Tm„ [eV] [prim.+N(z)+Pc(z)] 


ACT-A 


0.67 


0.61 


0.66 


0.61 


ACT-B 


0.67 


0.60 


0.65 


0.59 


ACT-C 


0.67 


0.51 


0.64 


0.49 


ACT-D 


0.67 


0.40 


0.61 


0.36 



TABLE V: Statistical uncertainty on total neutrino mass from cluster number counts and cluster correlation (similar to Table 
II) obtained from ACT-A, B, C and D SZ surveys (no LE). Shown are neutrino mass uncertainties (la) expected to be obtained 
from the respective planned surveys. We consider <tm„ assuming fiducial mass 0.1 eV. We show the results from CMB alone, 
CMB+cluster number counts N(z), CMB+cluster correlation Pc, and CMB+N(z)+Pc. Primordial CMB priors from PLANCK 
are considered. 
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a Mu [eV](LE) 


a Mv [eV][LE+N(z)] 


a Mv [eVl[LE+Pc(z)] 


a M ,[el / ][LE+N(z)+Pc(z)] 


ACT-A 


0.14 


0.14 


0.14 


0.14 


ACT-B 


0.14 


0.14 


0.14 


0.14 


ACT-C 


0.14 


0.14 


0.14 


0.13 


ACT-D 


0.14 


0.13 


0.13 


0.12 



TABLE VI: Same as in Table V but with LE priors. 



